Relaxation and Emittance Growth of a Thermal Charged-Particle Beam 



Tarcisio N. Teles0 Renato Pakter0 and Yan Leviijl 

Instituto de Fisica, UFRGS, Caixa Postal 15051, 
CEP 91501-970 Porto Alegre, Rio Grande do Sul, Brazil 

We present a theory which allows us to accurately calculate the distribution functions and the 
emittance growth of a thermal charged-particle beam after it relaxes to equilibrium. The theory 
can be used to obtain the fraction of particles which will evaporate from the beam to form a halo. 
The calculated emittance growth is found to be in excellent agreement with the simulations. 
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The understanding of physics involved in the transport 
of high-intensity charged-particle beams is of fundamen- 
tal importance in the development of a new generation of 
accelerators and electromagnetic wave generators to be 
used in applications such as heavy ion fusion, high-energy 
physics, communication, materials processing, and can- 
cer therapy. A very detrimental effect that may seriously 
influence the efhciency of such devices is a halo forma- 
tion and emittance growth of the beam [l|, 0, H, 0, 
These not only cause degradation of the beam quality, 
but may also be responsible for the activation of accel- 
erator channel wall and pulse shortening in microwave 
devices. Emittance growth is generally associated with 
the relaxation of initially non-stationary beam towards 
a more stable stationary configuration. The emittance 
growth can be calculated if the final stationary distribu- 
tion is known [1,[1,0|- However, the determination of this 
distribution is not an easy task [^, because parti- 
cles in an intense beam interact through long-range forces 
which prevent the system from relaxing to the true ther- 
modynamic equilibrium [l^, [HI, [13] ■ Instead these sys- 
tems get trapped in metastable states, lifetime of which 
diverges with the number of particles. To understand the 
properties of these states, one can not use the standard 
statistical mechanics and new non-equilibrium theories 
must be developed p^ . 

In this letter, we will present a theoretical framework 
which will allow us to accurately calculate the density 
and the velocity distributions of particles in the final 
stationary state achieved by a space-charge dominated 
beam focused by a uniform external magnetic field. Our 
approach is based on the theory of violent relaxation in 
gravitational systems modified so as to explicitly 

account for the effects of single particle resonances 
responsible for the halo formation [Tl| . The theory is ap- 
plicable to arbitrary initial conditions. In this letter we 
will show how the theory can be used to accurately calcu- 
late the density and the velocity distributions, as well as, 
to account for the emittance growth of a charged-particle 
beam launched with a thermal (Maxwell) velocity dis- 
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tribution. The predictions of the theory will be tested 
against the molecular dynamics simulations. 

The physical system considered here is an intense 
charged-particle beam of perveance K = Ni, / j^vlm 
— where c is the speed of light in vacuo, and q, m, and 
7h = [1 — (wz/c)^]^^/^ are the charge, mass, and the 
relativistic factor of the beam particles, respectively — 
propagating with an axial velocity VzSz through a mag- 
netic focusing channel enclosed by a cylindrical conduct- 
ing wall located at r = r^, [1, [l^ . The external focus- 
ing magnetic field is given by B = BqSz- It is conve- 
nient to work in the Larmor frame, which rotates with 
respect to the laboratory frame with angular velocity 
— <lBo/2"fi,Vzmc, normalized to Vz- In the Larmor 
frame, the external magnetic field produces a parabolic 
confining potential U{r) = Kzr'^/2, with the focusing 
field parameter Kz — Vl\/c^. The effective electromag- 
netic scalar potential between the particles ip, incorpo- 
rates both the self-electric and the self-magnetic fields, 
and B'*. This potential satisfies the Poisson equation 
with the boundary condition ipi'^'w) ~ 0, 

VV = ~(27rX/iVfc)nfc(r,s) (1) 

where, is the number of particles per unit axial length, 
r is the position vector in the transverse plane, and 
n(,(r, s) = Nb J /d^v is the transverse beam density pro- 
file, given in terms of the one particle distribution func- 
tion /(r,v;s). In the Larmor frame, the dynamics of 
the beam reduces to that of a two dimensional one com- 
ponent plasma with logarithmic interaction between the 
particles, confined by a parabolic potential U{r). The 
axial coordinate s = z = Vzt plays the role of time for 
this two dimensional system. 

We will suppose that the initial (transverse) distribu- 
tion of the beam is Gaussian in velocity space and is 
uniform in crossection, 

/o(r, v) = —^eir,^ - r)e-('^V2.^), (2) 

where is the initial mean square transverse velocity, 
and Tra is the beam radius. The quality of the beam 
is inversely proportional to the emittance, defined as 
= 4 < X >, for a stationary beam. For the 
distribution the emittance is Eq = 2crrm. 
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It will be convenient to discretize Eq. ^ into a p-level 
distribution 



(3) 



where p°(r, v) = Q{v — Vj^i)Q{vj — v)Q{rjn — r), and 
Vj and rjj are the maximum velocity and the amplitude 
of the level j, respectively, with vq — 0. For a perfect 
description of Eq. Q , an infinite number of levels {p — *■ 
oo) in ([3]) will be necessary. In practice, however, we 
find that a small number of levels is already sufficient 
to provide a very accurate approximation for the beam 
dynamics. For a given value of p, the optimal values of rjj 
and Vj can be obtained by minimizing the functional J- = 
J ifo — fpY'pYd'^'v, with the constraints on the kinetic 
energy and normalization. 
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where Ai and A2 are the two Lagrange multipliers. Min- 
imization of Eq. [4] yields the optimal parameters 
and {vj}. The many-body dynamics of systems with un- 
screened long range interaction is governed by the col- 
lisionless Boltzmann (Vlasov) equation. The distribu- 
tion functions which satisfy the Vlasov equation evolve 
in time as the density of an incompressible fluid. In par- 
ticular this means that the p hyper- volumes, 7(?7j) = 
/ ^(/0(r, v) - ?7j)d'^rd'^v, of the distribution Q will be 
preserved by the Vlasov flow [T^ . 

In reference [lj| it was argued that the stationary so- 
lution of the Vlasov equation could be obtained by max- 
imizing the coarse grained entropy, with the constraints 
imposed by the conservation of energy and the hyper- 
volumes of the p levels of the initial distribution func- 
tion. For matched beams and water-bag initial condi- 
tions, the resulting distribution was shown to be in excel- 
lent agreement with the molecular dynamics simulations 
[1^. However, for mismatched beams, the plasma os- 
cillations result in parametric resonances, which lead to 
a significant particle evaporation. After the relaxation 
process is complete, the stationary beam phase separates 
into a cold core, surrounded by a halo of highly ener- 
getic particles. For a water-bag initial condition [p — 1), 
it was shown that the core was very well described by 
a cold Fermi-Dirac distribution with the temperature 
T « TV/40, where Tp is the "Fermi temperature" of the 
beam. The halo was reasonably approximated by a step 
function, with energy range of one particle resonance. 
The full distribution function had the form of 



/(r,v) =/,(r,v)-KA(r,v) 
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For a p-level system, which is used to approximate the 
thermal distribution given by Eq. ([2]), a similar phase 
separation will occur. The form of the core distribution 



function can be obtained, once again, by maximizing the 
coarse grained entropy to yield 



with 



-/3j7je(r,v)+aj 



(6) 
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where the mean particle energy is e(r, v) = /2 + U (r) -t- 
^{r), and /? and {aj} are the Lagrange multipliers for the 
energy and the hyper-volumes conservation. The oscil- 
lations of the mismatched beam excite the parametric 
resonances resulting in a halo formation [l3| . The pa- 
rameter X determines the fraction of the particles which 
will evaporate to form the halo of the beam. The coarse- 
grained distribution can no longer preserve all the hyper- 
volumes of the original fine-grained distribution function, 
so that only the lower energy hyper-volumes will be con- 
served, while the particles from the higher energy states 
will evaporate to form a halo. We find that the halo can 
be modeled accurately by the distribution 

A(r,v) =xe(ec-e(r,v))^- 
xe(e(r, v) - ec)e(6H - e(r, v))e-^(^-^^) (8) 

The extent of the halo [l| is up to one particle resonance 
energy, e^. The low energy part of the halo distribu- 
tion is flat, while for energies e > = it decays 
exponentially with exponent 7^8. We can now, in prin- 
ciple, numerically solve Eqs.([Tl [5l [T]) and ([8]) to calculate 
the the stationary distribution function /(r, v) of the re- 
laxed beam. There is, however, one problem. Equations 
d?]) and ^ contain p-l- 2 parameter: /3, {aj}, and x- The 
conservation of energy, norm and of lower energy hyper- 
volumes gives us p + 1 additional equations. 



/ 



ydWv/(r,v) = l, 



(9) 



where 1 < j < p — 1 and eg is the average energy per 
particle of the initial thermal distribution. 
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There, however, still remains one missing condition nec- 
essary to uniquely determine the distribution function. 
For water-bag distributions, this condition was provided 
by the requirement that in the relaxed state, the core 
temperature is very low, T « Tp/AQ. It is difhcult, how- 
ever, to numerically implement this condition for p-level 
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distributions. On the other hand, if we discretize the 
original thermal distribution into only one level {p = 1), 
the condition T w Ti?/40 is easily implemented and al- 
lows us to uniquely close all the equations and calculate 
the relaxed distribution function [ISJ. We find, that al- 
though the core distribution is not well described hy ap — 
1 system, the halo part of the distribution is found to be 
quite accurate. This allows us to fix the value of x- Using 
this X, we can now improve the description of the core re- 
gion by including additional levels into the discretization 
procedure. To compare the predictions of the theory with 
the simulations, we calculated the number of particles in 
the interval [r, r -I- dr], N{r)dr = 2TTNi,rdr J d^v/(r,v), 
and the number of particles with velocities between v and 
V + dv, N{v)dv = 2nNbvdv J (i^r/(r, v), for various ini- 
tial conditions. The simulations are based on the Vlasov 
dynamics in which particles interact with the mean-field 
potential. This avoids the coUisional effects present in 
finite size systems, but which must vanish for one com- 
ponent plasmas in thermodynamic limit. The simulation 
code uses the Gauss law to calculate the mean electric 
field felt by each particle For axisymmetric beams 
studied in this work, this proves to be very efficient since 
the electric field at a radial coordinate r is determined 
simply by counting the total number of particles with 
coordinates smaller than r. Simulations were performed 
with 20000 particles. As can be seen from the Figs. 1 
and 2, the agreement between the theory and the sim- 
ulations is excellent. In the figures, the distances are 
measured in units of ^yeo/TtL and the velocities are in 
units of ^/eqUl- We have also defined a scaled perveance 
K* = K/fli^EQ and the mismatch parameter, fi = r,„/ro, 
which measures the deviation of the initial beam radius 
from the corresponding virial value tq — VK + Ao^ /Hl, 
for which the oscillations of the beam envelope are very 
small. In particular, we find that the discretization of 
the Gaussian by only 4- levels, already provides us with 
an almost perfect description of the core region. 
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FIG. 2: The density and velocity distributions: solid line 
is the prediction of the theory and points are the results of 
the molecular dynamics simulation. The scaled perveance is 
K* — 1 and mismatch is 50% (/i — 1.50). 



As a direct application of the theory developed above, 
we calculate the emittance growth of an originally ther- 
mal beam. This quantity is of fundamental importance 
for the design and develo pm ent of high intensity space- 
charge dominated beams [l6|. The calculations are per- 
formed for beams of varying scaled perveance K* and 
mismatch parameter /i. The results are compared with 
the molecular dynamics simulations. Once again, an ex- 
cellent agreement is found between the theory and the 
simulations (Fig. 3). 
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FIG. 1; The relaxed particle density of an initially thermal 
beam with scaled perveance K* — 1 and mismatch of 75% 
{fj, = 1.75). The points are the results of the simulations and 
the solid line is the prediction of the theory. Inset shows the 
exponential decay of the halo close to one particle resonance 
energy. 



FIG. 3: The final emittance of an initially thermal beam. 
The points are the result of molecular dynamics simulation 
and the fines are the predictions of the theory. 



To conclude, we have presented a theory which allows 
us to calculate the density and the velocity distributions 
of an initially thermal beam after it relaxes to the fi- 
nal stationary state. Comparing to the simulations, the 
theory is found to be extremely accurate, without any 
adjustable parameters. In particular, it can be used to 
calculate the emittance growth and the fraction of parti- 
cles which will evaporate as the beam evolves to its final 
stationary state. 
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